Isospin symmetry breaking and baryon-isospin correlations from 
Polyakov— Nambu— Jona-Lasinio model 
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We present a study of the 1+1 flavor system of strongly interacting matter in terms of the 
Polyakov— Nambu— Jona-Lasinio model. We find that though the small isospin symmetry breaking 
brought in through unequal light quark masses is too small to afl'ect the thermodynamics of the 
system in general, it may have significant effect in baryon-isospin correlations and have a measurable 
impact in heavy-ion collision experiments. 
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Signatures of phases of matter with deconfined color 
charges is under critical investigation for last few decades, 
both theoretically and experimentally. Quantum Chro- 
modynamics (QCD) is the formulation for first principle 
studies of strongly interacting matter. Along with the 
local color symmetry the quark sector has a few global 
symmetries. In the chiral limit for two light flavors u 
and d, we have global vector and axial vector symmetry 
S'?7v(2) (g) SUa{2). For non-zero quark masses, the axial 
symmetry S'[/a(2) is explicitly broken, while for non-zero 
quark mass difference vector (isospin) symmetry SUvi'^) 
is explicitly broken. At low energies the isospin sym- 
metry breaking (ISB) has relevance in many aspects of 
hadronic observables Apart from the quark mass 
difference, ISB effects may be brought in by electromag- 
netic contributions as well. Low energy tt — K scattering 
has been studied considering the inclusion of electromag- 
netic correction into the effective Lagrangian ISB of 
valence and sea quark distributions in protons and neu- 
trons has been studied in the chiral quark model 0, 0] ■ 
ISB may also have significant effect in the context of ex- 
istence of CP violating phase Some Lattice QCD 
investigation of the effect of unequal quark masses was 
done in Ref . Q and recently the effect of ISB on different 
hadronic observables were studied in Ref. 0, Q . Within 
the framework of chiral perturbation theory the isospin 
breaking effect in quark condensates has been studied 
considering niu ^ rrid and electromagnetic corrections as 
well, where the authors have given an analysis of scalar 
susceptibilities [ol. [lol|. 

In the context of high energy heavy ion collisions where 
strongly interacting matter is supposed to exist in a state 
of thermal and chemical equilibrium, the ISB effects have 
not been explored much. Fluctuations and correlations 
of conserved charges are important and sensitive probes 



for heavy ion physics. Most of the theoretical studies 
in this respect are in isospin symmetric limit (see e.g. 



11H27|). Here we present our first case study of ISB 



effect on fluctuations and correlations of strongly inter- 
acting matter within the framework of the Polyakov loop 
enhanced Nambu— Jona-Lasinio (PNJL) model. We dis- 
cuss the possible experimental manifestations of the ISB 
effects based on quite general considerations in the limit 
of small current quark masses. 

Formalism: In the last few years PNJL model has ap- 
peared in several forms and context to study the various 
aspects of phases of strongly interacting matter (see e.g. 
\WL [Tot). Here we use the form of the 2 flavor PNJL 
model with the thermodynamic potential studied by us 
in 2^, 211, but with a mass matrix, 



m = mil2x2 - m2T3 

mi — 7712 

mi + 771,2 



rriu 
ma 



where, l2x2 is the identity matrix in flavor space and T3 is 
the third Pauli matrix. While a non-zero mi breaks the 
chiral S'C/yi(2) symmetry explicitly a non-zero TO2 does 
the same for the isospin 5J7y(2) symmetry. For the case 
at hand m2 is also the difference of constituent masses of 
u and d quarks. 

We restrict ourselves to the mean field approximation 
and investigate the behavior of different charge suscep- 
tibilities. Once the pressure P is obtained, the n*'' or- 
der diagonal and off-diagonal susceptibilities are given by. 



and X 



XY 



with 



difix/TY' difix/TYdifiY/Ty 
i + j — n. These susceptibilities are related to the fluc- 
tuations and correlations of arbitrary conserved charges 
X and Y with corresponding chemical potentials fix and 



2 



IjLy- Given a two flavor system the global charge conser- 
vation is expected for baryon B, isospin / and electric Q 
charges. Since in the isospin symmetric limit the B — I 
correlations vanish exactly, we shall study this correla- 
tion with explicit isospin breaking. 

Results and Discussions: Here we consider the average 
quark mass mi = (m„ + md)/^ fixed at 0.0055 GeV and 
study the effect of ISB with three representative values 
of m2 = {md — m„)/2. The parameter set in the NJL 
sector has been determined separately for the different 
values of TO2. The differences in the parameter values 
were found to be practically insignificant. The bulk ther- 
modynamic properties of the system expressed through 
pressure, energy density, specific heat, speed of sound etc. 
did not show significant dependence on m2. Even the di- 
agonal susceptibilities were almost identical to those at 
the isospin symmetric limit. However, interesting differ- 
ences were observed for the off-diagonal susceptibilities 
in the B — I sector. We first discuss the fiB = case. We 
shall consider jij = in this work. It should be noted 
that even for /i/ = 0, a non-zero m2 would generate some 
non-zero isospin number. 
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FIG. 1. Second order off-diagonal susceptibility in B — I sector 
at hb = 0. 

In FigH] the second order off-diagonal susceptibility 
Xii , is plotted against T/Tc for different values of m2. 
Here Tc is the crossover temperature obtained from the 
inflection point of the scalar order parameters - the mean 
values of chiral condensate and Polyakov Loop 
As expected we find Xii = for ™2 = 0- For non-zero 
7712 the non-monotonic behavior observed here can be un- 
derstood as follows. At low temperatures the fermionic 
excitations are suppressed due to their large constituent 
masses as well as the suppression due to the Polyakov 
loop. On the other hand at high temperatures the mass 
difference 7712 become insignificant compared to the tem- 
perature scale. Therefore only at some intermediate tem- 
peratures one can expect a non-zero Xii ■ The peak value 
appears very close to Tc- 

The sensitivity of Xii on m2 is clearly visible. An 
exciting feature observed here is the almost linear scaling 
of Xii with m2 . This is shown in the inset of FiglTJ 

At the fourth order we have the off-diagonal susccpti- 
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FIG. 2. Fourth order off-diagonal susceptibilities. 



bilities, Xii ^ X31 ^-nd X22 ■ For fiB — the T dependence 
for the first two quantities and their linear scaling with 
7712 is shown in FigO For X22 '^^ found no such scaling 
behavior and the 7^2 dependence was insignificant. 

Qualitatively one can understand the behavior of X13 
and X31 by noticing that these are correlators between 
the fiuctuations xi of isospin and X2 baryon number re- 
spectively with the B — I correlator xfi ■ In our earlier 



studies |19l - l2l| we found that the both the fiuctuations 



increase monotonically with temperature. Here we found 
in Fig [1] above that with increase in T, Xii first increases 
and then decreases with a turning point close to T^- Thus 
Xii is correlated with X2 ^-^id xf below Tc and is anti- 
correlated above. Therefore one can expect that the cor- 
relations given by Xii ^^nd Xsi ^-r^ positive and negative 
respectively below and above T^. 

To understand the presence of 7712 scaling for some cor- 
relators and absence in others we first note that the dif- 
ferent B — I correlators may be expressed in terms of 
those in the fiavor space. The corresponding relation be- 
tween the chemical potentials are ^„ = \^lB + ^-nd 



This implies. 
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(1) 



The fiavor diagonal susceptibilities can be expanded in 
a Taylor series of the quark masses around to„ = m,d = 0. 
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(2) 



n=0 i=0 



where, 



1 



dral^dm-'. 



efficients, with i + j 



are the Taylor co- 

n and / £ u,d. Here Aqq and 
are respectively u and d flavor susceptibilities in the 
chiral limit; hence they are equal. Moreover, response of 
X2 to a change in m„ and that of X2 to a change in 
fTid {mu) are identical in the chiral limit. Thus we have 
= a^ j, V i, j. Therefore we get, 

X2 order) - Xjl'^*'' order) 

n 
i=0 



)• 



(3) 



where = a"„_j = a^_jj. It is clear that for any 
given n and i, the R.H.S. contains a factor {rud — rriu). 
Therefore Xii (EqUl) is proportional to m2 if the higher 
order terms are sub-dominant. This is what we observed 
for the range of 7712 considered here. 

For the higher order correlators one can similarly write. 











X22 = 





/Ltd\ 



i2i) 



•2X3^) 



(4) 
(5) 
(6) 



For all these quantities the first two terms on R.H.S. were 
found to be dominant. Considering again the Taylor ex- 
pansion in quark masses, X13 and X31 were found to be 
proportional to m2. No such proportionality was found 
for X^l 

Finally let us discuss Xii for non-zero /^s. In Fig[3]we 
present the results for three different temperatures - one 
above, one below and one close to the crossover temper- 
ature Tc- For large T, jiB effects are small and Xii while 
remaining positive slowly approaches zero with increas- 
ing jiB ■ Close to Tc Xii changes fast with increase of fiB 
and in fact changes sign and then slowly approaches zero. 
This effect becomes much more prominent at lower tem- 
peratures and at the first order phase transition bound- 
ary changes sign almost discontinuously. The negativity 
of the correlator indicates that the isospin number de- 
creases while increasing the baryon chemical potential. 
This is not expected for an ideal gas of quarks even with 
non-zero masses. Given m2 > 0, we have m„ < and 
therefore it is expected that Xii > 0- It is thus apparent 
that interactions play a major role in the high density 
matter. 

An amazing fact remains that the scaling of the corre- 
lators with m2 survives for all conditions of T and ^b- 
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FIG. 3. Second order off-diagonal susceptibilities along chem- 
ical potential for three different temperatures. 



This is shown in the insets of FigI31 Now expanding xfi 
in a Taylor series in hb about hb = we have, 



xf/(MB) = xf/(0) + ^X3^/(0) 



§xf/(o) 



,(7) 



In the above series odd order terms vanish identically as 
X^^ is CP even. Since Xii if^s) on the L.H.S scales with 
m2, the same can be expected to hold true individually 
for all the coefficients on the R.H.S up to any arbitrary 
order. Incidentally the first two Taylor coefficients have 
already been shown to respect this scaling in Fig[T] and 



Fig 2(b) respectively. 

Connection with Heavy Ion Collisions: Correlation be- 
tween conserved charges, is an experimentally measur- 
able quantity obtained from event-by-evcnt analysis in 



4 



heavy-ion collisions. The sign change of Xii high fiB 
may be an important signature of phase transition in this 
region. This will of course need fluctuations to freeze out 
close to the critical value of fiB- 
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FIG. 4. Ratio of B— I correlation to baryon number fluctua- 
tion at = 0. 

To compare with experiments it is often useful to con- 
sider ratios like R2 ~ Xii 1X2 ■ The temperature varia- 
tion of R2 is shown in FigH) It decreases monotonically 
to zero above Tc. Similar behavior can also be found 
with variation in /.t^ and the ratio goes to zero outside 
the phase boundary. A systematic study of this ratio can 
thus indicate how close one could approach the phase 
boundary in heavy-ion collisions. 

Also using the 7712 scaling one can estimate the mass 
asymmetry of constituent fermions in a physical system 



as, 7712 



cxpt 



i^^^ X 7712 P^JL. To the best of our 



knowledge this is the first theoretical formulation that 
indicates that quark mass asymmetry in thermodynamic 
equilibrium can be directly measured from heavy-ion ex- 
periments. There still remains the question whether the 
isospin asymmetry brought in through QED effects may 
disturb the scaling. 

Another important point to observe is that for frac- 
tional baryon number of the constituents Xii > Xsi ^ 
and vice versa for integral baryon number i.e. for pro- 
tons and neutrons. From Figl2] we sec that the former 
inequality persist well below Tc- This may well be an 
artifact of PNJL model. Thus it would be important to 
see corresponding results from Lattice QCD. Enhanced 
statistics of present and future experiments may make it 
possible to measure this extremely sensitive probe. The 
direction of the above inequality would be important in 
deciding if partonic matter may have been produced in 
the experiments. 
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